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ABSTRACT : 

The predicted dependence of sound transmission on the statistics 
of the randomly-r ough interface between dissimilar fluids has been 
studied by use of the Helmholtz Integral (J. Acoust. Soc. Amer . 51 f 
1083-1090 ( 1972 )). The predictions have been verified for radiation 
from a helicopter hovering, and slowly moving, over the sea , for fre- 
quencies to 1000 Hz for a wide range of surface acoustical roughnesses, 
R = k|cT^[ ( 02 / 01 ) cos 0^-cos ©2^^* roughness parameters a,k,c and 

0 are the rms height of the surface, propagation constant, speed of 
propagation and angle with the normal, respectively; subscript 1 re- 
fers to air and 2 to the water . 

When the surface is mirror -smooth , the transmitted sound pressure 
is due to the change in divergence and change in impedance at the 
interface. For low roughness, R < 1, the mean square transmitted 
pressure becomes decreasingly coherent and is a function of frequency; 
its magnitude is decreased by the factor e"^ compared to the perfectly 
smooth surface. For R ^ 1 the incoherent component dominates and the 
transmitted pressure depends also on the correlation length of the 
sur face displacements . 

The transmission change of sound pressure as a function of fre- 
quency is presented for several conditions of an SH3-D helicopter 
hovering and flying over or near an array of microphone and sonobuoy 
hydr ophones . 
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I. THEORY 



The Helmholtz Integral has been used, together with the Kirchhoff 

Assumption, to predict the dependence of the coherent (phase cancelled) 

component, and the incoherent component of airborne sound transmitted 

through a randomly rough ocean surface^. The solution is given 

simply in terms of the rms height. O’, and the correlation length, L, 

of the surface. The development is very similar to that used by many 
2 -7 

others to study the scatter ing of sound from a rough surface. 

Reference (1) describes various forms of the problem. We sum- 
marize these results by using square brackets [] to repeat the appro- 
priate equations of that work. We recall the cases: 

A) Infinite Plane Waves Incident on a Smooth Plane Interface in which 
the transmitted pressure ratio at the surface is, 

P 

^Is 



T = -21 = 1 * ~ ,1^1, 

Pio p^c^/cosQ^ + p^c^/cosQ^ 



[3] 



11 ' 



subscripts 1 refer to the medium of incidence. 

2 refer to the medium of transmission. 

0 = angle measured with normal 

p = density 

c = speed of sound 

= acoustic pressure, incident from medium 1, at the 
surface 



P2s 

R 

T 



= acoustic pressure, in medium 2, at the surface 
= pressure reflection ratio for plane waves 
= pressure transmission ratio for plane waves 



5 



B) Aperture-limited Plane Waves Incident on a Smooth Plane Interface 
(see Fig. 1 ) for which the Helmholtz Kirchhoff solution for a large 
aperature (AA K) is given by equation [6a]: 



[5] 



= p 



iK-r 



2m 






dxdy 



[ 6 ] 



where 



= p^ F sine K X sine K Y 
^2ra X y 



2m 



■i2p^0*^2® 






4TTr , 



AA cos 9 , 



AA = 4 xY 



1 k r * 



ID = 



~^ 2^2 -ik -r 
e 2 

e 



r = 



^2 ^2 

i^x + i^y + ^3^ » radius vector to the surface element 



r^ = distance from source to origin of ensonified area 

= distance from source to surface element 

r^ = distance from origin of ensonified area to receiver 

' ’ = distance from surface element to receiver 



sine K X = 

X 



s in K X 

^ 

K X 

X 



F = 1 + 



(C2/C1 sin - sin Q^) tan ©^ cos ©^ sin ©^ tan ©^ sin ©^ 



c^/c^ cos ©^ - cos ©2 



c^/c^ cos ©^- cos ©2 



^lA ^^^1^1"^^ 

P^Q = ~ — c y incident acoustic pressure at surface origin 



p^^ = acoustic pressure at unit distance from source 
p^ = acoustic pressure at the receiver 
dS = dxdy/cos y, area of a surface element 
C = C(x,y), surface height at (x,y) 
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GEOMETRY OF THE 
PROBLEM 




Figure 1. Geometry for Helmholtz -Kir chhoff theory for transmission 
of a plane wave through a randomly rough interface betwee 
two dissimilar media. 
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A 

n 



V = 
O' = 

= 

^2 = 

= 

= 

K = 
K = 

X 

K = 

y 

K = 



i^(-sin y cos 0')+ i 2 (-sin y sin 0')+ cos V, normal to 
surface element 

angle the normal makes with the z axis 

angle the projection of the normal onto the xy-plane makes 

with the X axis 
2TTf ^ ^ 

(i_ sin 0^+ i- cos 0^ ) , vector propagation constant of 

C^ 1 13 1 IT IT 

incident field 

2nf AAA 

(i sin 0 cos 0 + i^ sin 0^ sin 0«+ i cos 0^), vector 

^2 1 ^ 3^ ^ j3 ^ 

propagation constant of received field 

angle of incidence, measured from normal to interface 

angle that k^ makes with the z-axis 

angle that k^ makes with the xz-plane 

k, - k^ = 1 K + Lk + i k 
1 2 lx 2y 3z 

k, sin 0, = k^sin 0^ cos 0^ 

112 2 3 

k 2 sin 0^ sin 0^ 
k cos 0 - k cos 0 

X X ^ ^ 



K 



xy 



2 2 

(K + K ) 
' X y ' 



C) Aperture-limited Plane Waves Incident on a Rough Plane Interface 
(Gaussian PDF of heights, Gaussian correlation function) for which 
the Helmholtz-Kirchhoff solution is given by 

^P2P2 

(P2m) 

where R = ~ acoustical roughness for transmission 

(7 = rms height of surface 



= e"^ sinc^ K X sinc^ K Y + ^ 

X y AA 



„ " r,n 

e E r 

n = l nnl 



-K ^L^/4n 

Cio] 
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D) Radiation From a Point Source in Air Incident on a Smooth Plane 
Air-water Interface, for which the sound pressure transmission change 



for the mirror surface, TC^ (referred to is 



2m 



TC„ = 20 log^o — = - 20 



Is 



1 + f: 



D ®1 



-1 



H cosQ^ 



2cos9^ cosQ^ 



Cl4a] 



where D = depth = perpendicular distance from receiver to interface 
H = height = perpendicular distance from source to interface 



2ra 



= pressure at Snell angle 9 depth D, for a mirror -like surface, 



and for which the change referenced to the incident pressure at ground 
zero (P^2^ given by 



■^Siz = ^°®10 ®1 






E) Radiation From a Point Source in Air, Incident on a Randomly -rough 
Plane Air -water Interface. 

This last problem is, needless to say, the most difficult. The 
approximate approach derived in Reference 1 and used here (Section III- 
C-3) is similar to that in Reference 7. The ensonified area is divided 
into subareas; then all the coherent contributions and, separately, all 
the incoherent contributions, are added. The often incompatible criteria 
for each subarea are: 1) the incident pressure is approximately constant 
and the spatial variation of the incident phase is negligible (e.g. less 
than X/8) over the sub-area; 2) the linear dimension of the subarea is 
several surface correlation lengths in extent; 3) the correlation length 
is much greater than the acoustical wave length. The desirability of 
assumption (1) in adapting the plane wave solution is obvious. Assumpt- 
ion (2) states that the infinite integration that yields [lO] is accom- 
plished after "enough” correlation distances so that the spatial cor- 
relation has reached a "small" value. We will say more about this later 
(Section III-C-3). Assumption (3) is the requirement for our use of the 
Kirchhoff Assumption, in which the infinite plane wave impedance relation 

at the surface is assumed on a point-by-point basis. 
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For the point source we therefore start with Eq. [lo] for the 



transmission through an element of area AA which we now call , 

2 2 

multiply by (p ./p, .)> generalize by writing the sura for many 

such areas, (see Fig. 2). 






^2mi „2 



Is i 



" Plsi 



^2 2 
F. sine K -X. sine K .Y. + 

1 XI 1 yi 1 



2 2 TTL 

Po . F. -R. “ 

„ ^2mi 1 1 „ 

li — — : e 



1 P 



Is i 



A. 

1 



r" -K^ ,L^/4n 
1 xy 1 

e 



-n nn : 



n = l 



[ 15 ] 



We will assign i = 1 to the terras along the Snell (ray acoustic) 
direction and for simplicity we drop the subscript 1 for this most 
important direct ion . 

An approximate solution to the problem is obtained if we assume 
that there is no cross -coupling between coherent and incoherent com- 
ponents from each subarea. Then our problem becomes one of separately 
evaluating the terms 

<P2P2"^ , ^ ^P2S^ 



Is 



Is 



Is 



where the two terms on the rhs of Eq. [l6] represent the respective 
sums of coherent and incoherent terms on the rhs of Eq. [ 15 ]. 

The first term of Eq. [l5l is a sum of coherent components; the 
relative coherent transmission along the Snell direction is 






2s 



1 + 



Is 



Dc^cosO^ 

HCj^cosO^ 



(2 cosO^ cosO^) . 



[17] 



The incoherent term is obtained by summing the effects of the 
sub-areas. Adding this to the coherent term, and taking 10 log^^. 
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Source 




Geometry for point source above a rough surface. 

Figure 2, Geometry for Helmholtz-Kir chhoff theory applied to 
essentially plane sub-areas of a diverging wave as 
it passes through a randomly rough interface between 
two dissimilar media. 
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yields 
TC = 10 log 



10 



<P2sP2s*^ 



Is 



= 10 log 



10 



-R 



1 + 






^2 

Dc cos© \" 

Hc^cosS^ 1 COS»jCOS»2>' 



+ E 



TT A.L^ cos^©^ . F? S(R.K L) 
1 2rai 1 ' 1 xy ' 






C20] 



wher < 



;e the quantities and S(R^, K^^L) are functions of the angles at 

the area element, ©, . and ©_. (see Fig. 2 and Reference 1). 

li 2i ' ' 

Equation C20] states that, for very low frequencies (R « 1), the 

~R 

summation (incoherent) term will be very small, e "^1, and the 
transmission change will be given by the geometry and pc mismatch, 

(Eq. [ 14 ] • This pressure change may be positive or negative, i.e. 
a transmission gain or loss, respectively. As the frequency (and 
roughness) increases in the regime 0 < R < 1, the coherent transmission 
term will decrease but will still dominate the incoherent. For R 1 , 
the increasing incoherent term will begin to determine the magnitude 
of TC in Eq. [20] . This will cause an increase in transmitted pressure 
until a high enough frequency is reached (value depending on L/A^) for 
the S(R, K^^L) terra to decrease with fr equency^ and to again cause a 
greater transmission loss. 



II. EXPER IMENT 



A series of laboratory studies were undertaken as a prelude to an 
experiment using a SH3D helicopter as a noise source over the sea. The 
laboratory results and two graphs of the ocean data have been in- 
cluded in a paper titled ^’Spectral Characteristics of Sound Transmission 
Through the Rough Sea Surface”, by H. Medwin, R.A. Helbig and J.D. Hagy, 
Jr . , which has been submitted for publication in the Journal of the 
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Acoustical Society of America { es t iraated publ ica tion date Dec 1972). 

This technical report contains the complete presentation of the 
transmission results of the sea trials. 

The point source theory for a rough interface, Eq. C20], was tested 
by experiments at sea using the geometry as shown in Figs. 3 and 4. The 
”point^’ source was a Sikorsky SH-3D helicopter, at various heights, 
hovering or slowly moving on command over or near, an array of one 
microphone and several hydrophones, on command. The signal from the 
receivers was carried by cable to the nearby Floating Instrument 
Platform, FLIP, operated by the staff of the Scripps Institution of 
Oceanography. SQ^IPPS scientists made measurements of the sea surface 
height spectrum, concurrently with this experiment. 

Two series of experiments were conducted. FLIPEX I on 28, 29 May 
1970 and FLIPEX II on 17, 18 August 1970 both at a deep sea location, 

45 miles NW of San Diego. Only the results of FLIPEX II have been 
pr ocessed. 

A. Exper imental Details 

The SH3-D helicopter, viewed as a sound source, is characterized 
by a large number of strong harmonics of the very low frequency blade 
passage tone of the rotor, harmonics of the high speed counter-rotation 
tail propellor, and a broad band aerodynamic and machinery noise spec- 
trum, in which the harmonics are partially immersed. Only helicopter 
heights of 300 ft or greater are reported here; hovers at lower levels 
disturbed the sea surface and precluded comparison with theory. 

The hydrophones were parts of modified AN/SSQ-57 sonobuoys. The 
two modifications consisted of: a) disabling the RF modulator, which 
is normally used to radio the signal from the buoy; and b) changing the 
length of the hydrophone support so that it hung at depth 0 (or 10 ft). 
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Figure 3- Geometry of transmission into 
source in air. The Snell ray 
(Not to scale.) 



water from point 
is indicated. 
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FLIPEX I 





FLIPEX H 



CLUSTER I 




40 '^ 



CLUSTER H 







ZOO'6 



200' 6 



SIDE VIEW OF EXPERIMENT 
ARRANGEMENT. 



Fiaure 4. Experimental arrangement for ocean experiment. Noise 
Figu Trs 

or Cluster II. For FLIPEX II, reported here, the clusters 
were separated by 150 ft. Ocean and acoustic ^ta were 
recorded at Scripps Institutions Research Platform FL , 
300 ft upwind from Cluster I. 
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20 , 40 , 95 , 200 or 300 ft below the surface. The final configuration 
of microphone was simply a hydrophone placed on the bottom of a 4 ft 
diameter plastic Kiddy Pool. The colorful Kiddy Pool served not only 
as the floating holder for the microphone but also as an easily identi- 
fiable reference position by which to guide the helicopter . The modi- 
fied sonobuoy consisted of hydrophone, pre-amplifier, audio amplifier 
and transformer powered by a salt-water battery. Three sonobuoys were 
held together in a cluster. There were two clusters as shown in Fig. 4 , 
one 300 ft from FLIP, the other either adjacent to or separated from 
the first by an additional in line, 150 ft. The microphone and hydro- 
phones were hard-wired to a tape recorder in the instrument room on 
FLIP. The electrical cable was floated on one inch polypropylene line. 
Additional details of the mechanical, electrical and acoustical design 
and tests will be found in Ref. 11 . 

The signal from the sonobuoys was FM recorded on an 8 channel 
PI 6208 magnetic tape recorder, on board FLIP. One channel was re- 
served for voice annotation. 

The data recorded during FLIPEX II consisted of calibrated, time- 
varying, voltages from hydrophones at 5 depths and a surface microphone 
during helicopter hovers at heights 40, 600 and 750 ft and 5 knot flybys 
and fly-overs at 150, 300, 600 and 750 ft. 

Under the direction of Dr. Russ Davis of Scripps, time series wave 
height variations were recorded during both FLIPEX I and FLIPEX II. 
Whenever possible, such data were taken simultaneously with helicopter 
noise data. During FLIPEX II, four wave probes were used. Surface 
wave frequency spectra were generated by transforming the wave height 
variation time records for each device. Then mean square wave heights 
for several wave period ranges were obtained by integration of the 
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frequency spectra over appropriate frequency ranges. This processing 
was carried out at Scripps. 

B. Data Analysis 

1. Analogue Analysis 

The frequency range 100 Hz to 1000 Hz corresponded to an 
acoustical roughness range, R, of approximately 0.05 to 4.0 for nor- 
mally incident sound and for the typica 1 sea surface encountered 
during the experiment. Sixteen frequencies in this frequency range 
were selected for investigating the effect of a broad range of surface 
roughness on transmission and for comparison with theory. 

The analysis utilized a Hewlett Packard Model 3590A wave 
analyzer with a General Radio Type II 63 A frequency synthesizer used as 
an external local oscillator for the wave analyzer. The logarithmic 
output of the wave analyzer was used as the vertical input to a Varian 
Model F-lOO X-Y recorder. This arrangement was particularly useful for 
analyzing and displaying the helicopter passes since the resulting plot 
showed the ambient noise levels before and after each pass as well as 
the total noise during each pass. A series of plots made for several 
frequencies using a typical helicopter pass appears in Fig. 5. The 
large time variation in rms level even with the slow response time of 
the wave analyzer is evident. Signal voltage levels for the helicopter 
passes were determined by taking an eyeball average of the time vari- 
ations across the peak of the pass. 

Ambient noise voltage levels between passes were determined in the 
same way. Since the plots represented rms levels in dB, this eyeball 
average of dB levels was something less than the true average rms voltage 
level. The error here was on the order of 0.5 dB. Since the type of 
time variation was about the same for the surface microphone as it was 
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Figur e 



! 





. Strip chart records of sixty second samples of 
noise received in six 10 Hz bands as helicopter 
passed over the hO ft deep hydrophone. 
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for the submerged hydrophones, this error canceled itself out when 
taking the difference to compute transmission change. To obtain heli- 
copter signal voltage levels, the ambient noise voltage levels were 
subtracted off by using the white noise assumption. Signal sound 
pressure levels at the hydrophone were then determined by using the 
calibrated sensitivities of the transducers. 

The Kirchhoff assumption was made in order to convert from ob- 
served microphone SPL to incident sound SPL. That is, it was assumed 
that pressure doubling had taken place so that the incident SPL was 
obtained by subtracting 6 dB from the surface SPL. 

2 . Digital Analysis 

The major data reduction program for the transmission study 
was performed digitally. 

Figure 6 is a general self-explanatory flow chart of the 
analysis. In addition to the flow chart, the following statements will 
help to characterize the analysis: 

The sampling frequency was 2000/sec. Four channels, three 
hydrophones and one microphone, were sampled simultaneously. For each 
flyby or flyover experiment, data were digitized for peak sound levels 
when the helicopter flew directly over the mic and/or hydrophones. 
Additional data were taken in order to study changing sound levels 
within a time span which began approximately 50 sec before the helo 
approached the overhead position and ended approximately 50 secs after 
the overhead position. 

Ambient noise levels, digitized for corresponding hydrophone 
depths, were taken at least one half hour preceding or succeeding each 
experiment. Most ambient signals were taken within minutes of an 
experiment and spanned a time period of approximately 15 seconds. 
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After proper formatting of the tape, processing of the data 

12 . . 

was begun using a system of time series programs originally developed 
by the University of British Columbia. The first program in the series, 
program FTOR , employs the FAST FOURIER TRANSFORM to compute Fourier 
coefficients from selected channels of the digital tape. 

Program FTOR was run using varying block lengths of data in 
order to arrive at an adequate frequency resolution or bandwidth. After 
some experimentation using block lengths ranging from 2 secs to 1/8 sec, 
the decision was made to use blocks of 1024 samples (a duration of ap- 
proximately 1/2 sec), giving a frequency resolution of 1.95 Hz. 

Contiguous blocks of 2^^ data points were transformed for each 
experiment and written on a ’Courier Coefficient Tape^^. This tape was 
then read by program FCPLOT of the time series programs. The FCPLOT 
program provides the capability of averaging over only selected blocks, 
e.g., over blocks spanning 10 second periods, or over all blocks as 
was the case for ambient data . This program was revised to write power 
spectral densities for given time spans on magnetic tape for further 
analyses . 

From this tape of power spectral densities created by program 
FCPLOT, blocks of data corresponding in time were further averaged for 
different runs of a given experiment. For example, peak data were 
averaged for several runs and several ambients were averaged. Data 
from different runs were weighted according to the time span they re- 
presented. The averaged power spectral densities, along with their 
associated ambient levels, were read by program TC, a program which 
calculates sound pressure levels and sound transmission change and plots 
them against frequency and roughness, respectively. 

Ambient noise was subtracted from helo-plus -ambient levels. 

Points were labelled ^^bad'^ if (helo-plus-ambient) -(ambient 0, or if any 
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SPL differed by 20 dB or more compared to the SPL at a contiguous 
frequency . 

An average TC was calculated for every 9 contiguous frequencies 
(band of 9*75 Hz). This average was considered ^’good” if 5 points out 
of the 9 were considered ’’good’*. 

Because of many ’’bad” points caused by high ambient levels for 
f < 100 Hz, the ’’low frequency SPL” was calculated for the frequency 
band 100 to 200 Hz; the ’’overall SPL” was calculated for the frequency 
band 100 to 1000 Hz. 

C. Comparisons of Experiments with Theory 

During the course of the acoustical study, ocean height data were 
taken simultaneously at four probes, at a sample rate of eight times 
per second. The motions of FLIP were compensated for. The data for 
each probe were divided into four minute blocks, the spectra were 
computed, and then the spectra were averaged across these blocks. The 
wave height uncertainty for each of the runs is probably less than 5 %» 
The wave height frequency spectrum for each probe was integrated to 
yield the mean squared height. The values of the mean squared height 
at the four probes were averaged and this average was used for the 
theoretical calculation of sound transmission changes. A calculation 
such as described was carried out for the ocean wave system twice on 
17 August and twice on l8 August to yield values of rms heights. O’: 

15*5 cm, l4*9 cm, 13*6 cm and l4*5 cm. 

The general comparison with theoretical predictions is comprised 
of several single parametrical checks: 

1 . Variation of Transmission with Depth/Heiqht Ratio for Normal 

Incidence at a Very Smooth Surface R « 1 

When the surface is relatively smooth, R « 1, Eq. [20j tells 

us that only the coherent component of the transmission change is 
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significant. Further, if the surface is \ery smooth, R « 1, TC 
approaches (Eq* l 4 a ) and the transmission change will be due only 

to the geometry and the pc mismatch. The test is accomplished for the 
ocean experiment by calculating the average experimental TC for the 
low frequency band 100-200 H 2 . For this band the surface acoustical 
roughness for normal incidence sound (fly-over or hover) was R =0.15 

at 200 Hz and one-quarter that value at 100 Hz. The decrease in co- 

-R 

herent transmitted signal due to surface roughness, which goes as e , 
would therefore be expected to be at most 0.68 dB at the high fre- 
quency end of the band and O .16 dB at the low frequency end (100 Hz). 
This average expected loss of approximately 0.4 dB would not be dis- 
tinguishable in this experiment. On the basis of the above argument, 
all observations of transmission change for the low frequency band, 
100-200 Hz, are identified as occurring for R « 1 and are assumed to 
be equivalent to the values that would have been obtained for a mirror- 
like surface. 

Figure 7 shows the transmission change at normal incidence pre- 
dicted by mirror surface theory (R « 1) compared with experimental data 
taken in hovers and flyovers, at a height of 600 ft, for a decade of 
ratios of hydrophone depth, D, to source height, H. Only the low fre- 
quency band of sound frequencies have been used for this test. A 
positive value of TC represents an SPL at the hydrophone that is greater 
than the incident SPL at the surface. The average disagreement is 
approximately 2 dB. This small discrepancy could have been due to one 
or more of the following: the helicopter was not a proper point source: 
there was a difference in the absolute calibration of microphone and 
hydrophones; the pressure was somewhat less than doubled at the slightly 
rough surface. 
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2 . Variation of Transmission with Angle of Incidence at a Very 
Smooth Surface (R « 1) for a Constant Depth/Height Ratio 
Figure 8 presents low frequency data taken during fly-overs. 
Each datum point represents average over eight or l 6 blocks (each 
block is ^ second long) for cases in which the aircraft was either 
approaching (at 5 knots) or receeding. The geometry and the defi- 
nition of ^20 shown in Fig. 3- The assumption has been made 

(as above) that low frequency data for which (R « 1 ) should be com- 
pared with theoretical curves for the mirror surface TC,. The com- 

M 

parison is again favorable for the two D/H ratios, O.O 67 and 0.5? fot 
which data were available. The theoretical curve for D/H = 10, is 
drawn for interest. 

3 • Variation of Transmission with Sound Frequency and Surface 
Acoustical Roughness; Source over Receiver 

As the sound frequency increases, Eq. [ 20 ] predicts that the 
coherent term will dominate briefly and then for R > 1 the incoherent 
term will control the total transmission change. In order to calculate 
the expected transmission, Eq. [ 20 ] requires knowledge of the surface 
height correlation length L. This can be estimated from the known wind 
speed, by using the results of Cox and Munk to calculate first the rms 
slope and then using the assumption of a Gaussian Correlation Length, 
L, to calculate its value; this simple calculation yields the approxi- 
mate correlation length, L 150 cm for the ocean during FLIPEX II. 

A second decision must be made before Eq. C 20 ] can be used in the 
case of a point source; one must settle on the appropriate sub-areas 
AA^. In principle they should be chosen so that the incident sound 
wave has uniform phase and amplitude over the subareas, and the cor- 
relation has dropped to effectively zero at the boundary of the sub- 

area. Figure 9 shows the result of a computer study of the effect 
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of a) changing the correlation lengths, L, and b) changing the width 
of the circular rings AL which define the sub-area. The calculation 
is for a source at 600 ft and a receiver at 20 ft for a sea of rras 
height, Cr = 13 . cm and correlation length, L ^ 150. cm (= 5 ft). It 
is observed that only the higher frequencies (R > 1) are controlled by 
these changes (that is, only the incoherent component is modified). 
Changing the correlation length from L = 2.5 ft to 5 ft to 10 ft with 
AL = 3b has the major effect on the transmission; changing the size of 
the sub-area rings from AL = L to AL = 3L or AL = 5L, for the appro- 
priate L = 5 ft has a smaller effect on the predicted transmission change, 
TC. Theoretical calculations of TC to be presented, below, were based 
on L = 5 ft and AL = 5L which fulfills assumption (2) of Section I-E* 

It is not possible to simultaneously fulfill assumptions (1) 
and (3) of Section I-E for all the wave lengths studied in this exper- 
iment. Figure 10 presents the experimental data, analyzed by digital 
computer (solid line), and analogue playback, (circles). The two 
methods of analysis are in tolerable agreement. The low frequency ex- 
perimental values of TC greater than +6 dB are clearly due to ambient 
noise underwater. The theoretical predictions (dashed line) appear to 
have approximately the correct dependence on roughness (or frequency 
squared). However theory appears to underestimate the amount of sound 
transmitted by roughly 5 dB for both the hydrophone at 20 ft depth and 
the one at 300 ft. Perhaps 2 to 3 dB of this discrepancy is the same 
as observed for R « 1 (Figure 7) and the explanation of section (1) 
could be repeated. However, it is probable also that the simple theory 
for adding sub-areas that was developed in reference 1 underestimates 
the contributions of the incoherent sound in the region R ^ 1. 
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The transmission change for normal incidence sound is seen to 

be only slightly dependent on acoustical roughness (frequency squared). 

For example the experimental curves could be crudely summarized as 

showing a constant TC ^ 2 dB. It is this small spectral variation for 

near -overhead sources that has permitted some researchers to ignore 

the frequency dependence, completely, and to assume that the mirror 

value TC^ Eq. [l4] as valid for all frequencies. This fortuitous 

situation is due, in fact, to large incoherent contributions from non- 

Snell directions which compensate for the decrease in transmission 

of the coherent component, as the frequency is increased. 

The jaggedness in the experimental data is, of course, due to 

inadequate averaging time being available. Much of the data represent 

averaging over 10 to 20 seconds, whereas 10 to 20 minutes would be 

necessary to reduce adequately the statistical variations. 

Additional experimental curves of transmission change for the 

helicopter over the hydrophone are given in Figures 11-25. Figures 

11-13 represent hovers nominally over the hydrophones at 20, 40 and 

300 ft respectively. In fact, the helicopter pilot had difficulty 

seeing the float marker and maintaining position during the hover at 

600 ft altitude and unintentional off-sets of as much as 100 ft, may 

have been present. Our study of the playback gave evidence of the 

helicopter drift and permitted the extraction of the 24 second analysis 

presented in Figures 11-13 which we believe represents the true over - 

-R 

head position. The coherent, e , dependence for 0.2 < R 0.5 is fairly 

clear in these cases because of longer averaging time; the straight 

-R 

line starting at TC = 0 shows the predicted slope of an e dependence. 
The data for R < 0.2 are unreliable because of strong and variable 

ambient noise levels. The hover analysis yields virtually the same 
TC data as the fly-over analysis shown in Figure 10. 
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Figures l4 and 15 represent 5 sec and 3 sec analyses for source 
at 600 ft over the hydrophone at 20 ft. The 150 ft offset of the 
source relative to the microphone would have a small effect on the 
incident levels; an inverse-square correction has been made for this. 

The two curves show the type of variation of TC that can take place 
due to a different sea surface at the time of analysis and due to an 
inadequate time of analysis. Some of the jaggedness is eliminated 
when the 8 sec average is taken; this average was presented in Figure 10. 
Figure l6 represents an 8 second average for the overhead source at 
600 ft and receiver at kO ft; it should be compared with the curves 
on Figure 10 from which it was omitted to avoid crowding. 

Figure 17 shows the transmission change for a pass directly 
over Cluster II containing microphone and 95 ft hydrophone. It too 
can be compared with curves on Figures 10, from which it was excluded 
to avoid crowding. 

4 . Variation of Transmission with Sound Frequency and Surface 
Acoustical Roughness; Source Offset From Receiver 
When the source is not directly over the receiving hydrophone, 
our theory again gives a fairly good prediction of experiment. Figure 
l8 shows the comparison of theory and experiment for an offset of 
approximately 250 ft, resulting in an angle from hydrophone, 

85 ^ for the 20 ft receiver and 40^ for the 300 ft deep receiver. In 
the case of the 20 ft hydrophone, simple mirror theory would predict a 
transmission loss of more than 15 dB for all frequencies. The theory 
of Eq. C20] comes closer to the mark because it accounts for the in- 
coherent contributions which become important at R > 1. 

The region R < 0.5 is of particular interest because of the so- 
called ’’evanescent”, ’’inhomogeneous”, or ’’lateral” wave that 

13-14 

theorists predict would be refracted through a smooth interface. 
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TC tDB) 




EKP^ LI - HOUER 

HlLO RLTJTUDt: 6 00. 0 FT- HVORDPHCNE DEPTH 20-0 FT. 

HELD OFFSET FRDH HJC FLORT - URRJED 
HELD OFFSET FRDH HVQRO FLDRT - HRRIED 
THETRI -0.00 DEC. THETR2 --0.00 DEC. 

HO. OF ELKS - 46 I i^2 SET ELKS) 

Figure 11, Transmission change in SPL as a function of surface acoustical roughness 
for helicopter hover at 600 ft over hydrophone at 20 ft depth. 
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EKP. ii - HQUER 



HELO RLTlTUDt; -■ GOO. 0 FT- HVDROPHOHF DEPTH -- 40- P FT. 

HELD DFF5ET FROM H]D FLaPT - URRIED 
HELD OFFSET FROM HVDRG FLORT - URRIED 
THETRl ---0.00 DEG. THETR2 --0.Q0 DEG. 

HO. OF BLK5 -- 40 11^2 SEE BLkSI 

Figure 12. Same as Figure 11 but with hydrophone at 40 ft depth. 
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EKP. il - HOUER 



HELD RLTITUDE nOO-O FT. H^DROPHGHE DEPTH - 

HELD OFFSET FROM M]0 FLORT - UPRIED 
HELD OFFSET FROM H^DRO FLORT - URRJED 
THETRl ---0.00 DEG. THETR'2 ---0.00 DEG. 

MO. OF BLK5 - 40 iU2 SEE BlKSJ 

Figure 13. Same as Figure 11 but with hydrophone at 300 ft depth. 
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TC ^0B^ 




17 - RUN 16 - FL^BN 

HL'La RLTJTUDE --- 600-0 FT. HVORaPHONF DEPTH - 20-0 -T. 

HELD DFFSET FROM HJC FLQRT - IFO FT 
HELD DFFSET FRDH H'-'ORDU) - 0 FT- 
THETRl -0.00 DEG. THETR2 - O-Ofl DEG. 

MD. GF BlKS ■■ iO '1/2 SEG BLK5) 

Figure 14. Transmission change as a function of surface acoustical roughness for a 
helicopter fly-by at 600 ft over hydrophone at 20 ft depth. 




.17 - RUN 19 - PLVRN 

HP.lC RL.nTuDL ■- bOD- D PT. H^DRGPhnME DEPTH - 2 C -0 

HEia OFFSET FROM HID FLORT 150 ^'T 
HELO OFFSET FROH H'-'QROU) - 0 FT. 

THETRi -U.D;] DEO. theTR'2 0. 00 DEO. 

MG. OF BLK5 ' 0 1I"2 SEE DLK5) 

Figure 15. Same as Figure l4 but a different fly-by. 
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EKP. L? - FLVBE - RUG. OUER 2 RUNS 



HELD RLTITUDE -- 600.0 FT. H^nHOPHOHE DEPTH - 40.0 

HELD OFFSET FROM KUE FLORT iSD FT 
HELO OFFSET FROM HF-DROUJ - 0 FT. 

THETRL -O.QO OEG. THETR2 --'0.00 DEG. 

NO. OF DLKS -- 16 [ 1^2 SEE BLK5J 

Figure 16. Average transmission change as a function of surface acoustical roughness 
for two helicopter fly-bys at 600 ft over hydrophone at 40 ft depth. 
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13 - FL^nUtR - PERK PLUS 0 SECS 

Hf:ua RLTiTUDL - BOO-O FT. H^QRQPHONE DEPTH -- 35. D FT. 

HELD OFFSET FROH H!C FLORT - 0 FT. 

HELD OFFSET FROH H^DROdl) --.0 FT. 

Tf-MR: .^,uu QEG. T^^rTR^ -ij.Dj Rf'E;. 

HO- OF BLKS - 22 lid SEC BLKSI 

Figure 17. Transmission change as a function of surface acoustical roughness for 
helicopter fly-by at 600 ft over hydrophone at 95 ft depth. 
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Our theory has not considered this component. The fact that the low 
frequency sound levels underwater are higher than predicted by our 
theory (see also Fig. 8 with D/H 0,067) may be evidence that a 
later ial wave exists for rough surfaces, too; unfortunately, it may 
also be evidence of the strong low frequency ambient noise that was 
present during the experiment and/or refraction to the 300 ft hydro- 
phone due to a strong (9C^) negative temperature gradient that existed 
between 60 ft and 240 ft. 

Figures 19-21 show transmission changes as a function of 
roughness for a helicopter pass directly over the microphone but 
offset by 150 ft from the string of three hydrophones of Cluster I at 
depths 20, 40 , 300 ft. The helicopter was moving perpendicular to 
the line of Clusters I and II, when these offset data were collected. 

The next series of graphs. Figures 22-25, show spectra of TC 
at various times during 5 knot fly-overs at 600 ft altitude, using 
hydrophones at depths 20, 40 and 300 ft. Now the helicopter moves 
perpendicular to the line of Clusters I and II so as to pass over 

the hydrophones. In Figure 22 we see three curves for the 20 ft 

hydrophone representing (from bottom to top): source 50 sec (4l0 ft) 
offset and approaching; source 30 secs (250 ft) offset and approaching; 
source overhead. 

Figures 23 shows the spectra of TC for the receding helicopter 
as it continues to fly slowly along the path perpendicular to the line 
of the microphone-hydrophone clusters: passing over the 20 ft hydrophone 
in the top curve; 40 sec (330 ft) past the hydrophone cluster in the 

middle curve; and 50 sec (4l0 ft) past in the bottom curve. 
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13 - FLVaUER - PERK PLUS 0 SECS 



HllQ RLTJTUDL -- 600.0 FT. H^CROPHOHE DEPTH -- 

HELD DFESET FROh HiC FLORt 0 FT. 

HE^O OFFSET FROH ^^^DROdl -- 150 FT- 
THETPl --ll 13 DEO. THETP2 ■•57.09 DEG 

NO. OF BLK3 ■ 22 dd SEG BLKS) 



20.0 FT. 



Figure 19. Transmission change as a function of surface acoustical roughness for 
helicopter offset by 150 ft from 20 ft deep hydrophone. See text for 
details. i 
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EKP. [3 - PL^aUER 

Hf'LQ PLTJTUCt -■■ 600-0 FT. H>^0RaPHC]NF DFPTH - 40-0 -T. 

HlLQ CFFSFT from MIC. p-LOPT 0 FT- 
PELO OFFSET FROF H'^OROUJ - ISO FT- 
THETPT -g.OJ QFC. THETP2 ---'40-02 GEG. 

HQ. QF BLK5 23 ' SEC RLK5I 

Figure 20. Same as Figure 19 but with hydrophone at 40 ft depth. See text for details. 
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13 - FL^aUER 

HLLQ PLTJTUCE - 600- G FT. Hi/npnpHaH'L DL'Pl'H -■ jGC.r '■ 

HFLQ GFFSFT FRn[^ ^^IC, FLGPT ■- 0 FT. 

HFLQ QFF5TT FRdM ^^^nRDUJ - 150 FT. 

THFTRl -4.3b DLT;.. THETP2 -IG-iP PEb. 

RQ. OF BL.K5 23 i 1-2 SEE RLKSI 

Figure 21. Same as Figure 19 but with hydrophone at 300 ft depth. See text for details. 
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Figure 22. Transmission change as a function of surface acoustical 

roughness for H = 600 ft, D = 20 ft at different positions 
of flyby. See text for details. 
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Figure 23 . Same as Figure 22; see text for details. 
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Figure 24 shows TC under the same conditions as described in 
the previous figure except that hydrophone is 40 ft below surface. 

The top curve is for the overhead helicopter; the middle curve for 
the source receding and at I 65 ft from overhead; the bottom 

curve for the source receding and at 250 ft from overhead. 

Figure 25 presents TC for the same conditions of flyby except 
that the hydrophone is 3 OO ft below the surface. The top curve is 
for the helicopter overhead position; the middle curve is for the 
receding helicopter now at a distance of 250 ft from overhead; the 
bottom curve is for the receding helicopter at 4l0 ft from overhead. 

5 . Other Parameters 

Temperature refraction was an unavoidable variable during the 
experiments reported here. Figure 26 shows the BT plot typical during 
FLIPEX II. It is believed that the sound pressures at the 20 ft 
hydrophone were unaffected by the existing gradient. Nor would the 
results at normal incidence be affected for the hydrophones at any 
depth. On the other hand the TC curves at non-normal incidence for 
hydrophones below 20 ft are only approximate because we have not 
corrected the data for the temperature refraction effect. 

III. CONCLUSIONS 

The Helmholtz-Kir chhoff theory, applied to transmission from air 
into water, provides useful predictions of the dependence of trans- 
mitted sound pressure level on frequency, surface roughness and source 
and receiver positions. For low roughness, R < 1, transmission is 
essentially the same as for a smooth surface, but decreased by the 
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TEMPERATURE (®F) 




BT DATA 

I72I30T AUG. 1970 

Figure 26. Bathythermograph at time of experiments. 
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factor e . For larger acoustical roughnesses, R ^ 1, it is possible 
to make predictions that come close to reality provided the ensonified 
surface is judiciously divided into sub-areas • The experimental data 
are not able to clearly demonstrate the presence of an important 
lateral wave component for low roughness surfaces. 

Greater accuracy for the rough surface transmission problem will 
require a theoretical approach that takes into account the lateral 
wave in the case of a rough surface and which considers interactions 
between sub-areas or deals with the affects of diverging waves from 
the beginning of the development. 
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The predicted dependence of sound transmission on the statistics of the 
randomly-r ough interface between dissimilar fluids has been studied by use of 
the Helmholtz Integral (J. Acoust. Soc. Amer . IO 83 -IO 9 O (1972). The pre- 

dictions have been verified for radiation from a helicopter hovering, and 
slowly moving, over the sea, for frequencies to 1000 Hz for a wide range of 
surface acoustical roughnesses, R = ( C 2 /c]^)cos 9^-cos © 2 ^^* rough- 

ness parameters CT,k,c and 9 are the rras height of the surface, propagation 
constant, speed of propagation and angle with the normal, respectively; sub- 
script 1 refers to air and 2 to the water. 

When the surface is mirror -smooth, the transmitted sound pressure is due to 
the change in divergence and change in impedance at the interface. For low 
roughness, R < 1, the mean square transmitted pressure becomes decreasingly 
coherent and is a function of frequency; its magnitude is decreased by the 
factor e“^ compared to the perfectly smooth surface. For R ^ 1 the incoheren 
component dominates and the transmitted pressure depends also on the cor- 
relation length of the surface displacements. 

The transmission change of sound pressure as a function of frequency is 
presented for several conditions of an SH3-D helicopter hovering and flying 
over or near an array of microphone and sonobuoy hydrophones. 
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